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^ ■ This is an extended summary of the two parallel sessions held at MG11: PPN1 "Strong 

■^J- ' Gravity and Binaries" (chaired by L.B. and L.G.) and PPN2 "Post-Newtonian Dynamics 

£f*) ' in Binary Objects" (chaired by G.S.). The aims and contents of these sessions were 

f *^ , close to each other and overlapping. It is natural to review both sessions in one joint 

i contribution to the MG11 Proceedings. The summary places the delivered talks in a 

f**^ ' broader perspective of current studies in this area. One can find more details in individual 

' contributions of the respective authors. 

o ' 

1. Introduction and overview 

| | The current strong interest toward binary compact objects is largely driven by the 

imminent observation of gravitational waves by the presently operating or, more 
^ , likely, by the upcoming advanced detectors. Many astronomical systems can emit 

' some amounts of gravitational waves. However, it follows from the most general 

j_j ■ theoretical considerations that the amount of gravitational radiation is maximized 

when two massive chunks of matter are moving with respect to each other with 
relativistic speeds. In the Cosmos, this situation occurs naturally in tight binary 
systems involving the most compact objects presently known — neutron stars and 
black holes. 

This explains the special attention to the massive compact objects orbiting each 
other at the late stages of their relativistic evolution — the so-called inspiralling 
compact binaries. 1 ' 2 It is argued 3 ' 4 that the inspiralling pairs of stellar mass black 
holes will probably be the first sources directly detected by the ground-based laser 
interferometers. The theoretical description of such binaries is usually being based 
on the Post-Newtonian (PN) approximations to general relativity. 5-7 

In PN gravitational-wave studies it is often sufficient to treat black hole as a 
highly compact massive object, independently of whether the object possesses the 
general-relativistic event horizon or not. In the Newtonian limit, when formally 
c — » oo, the Schwarzschild radius r g — 2Gm/c 2 shrinks to zero, so the black hole 
can be viewed as a point- like particle endowed with a mass (and possibly a spin). 



1 
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The PN approximation is expected to provide a good description of black holes 
and their relative motion in an expansion series with c — > oo, that is, when the 
dimcnsionless ratios, such as r g /r and v/c, are sufficiently small. 

The masses of black holes are assumed to be ranging from stellar masses (up to 
a few tens of M©) to supermassive black holes (SMBH) in nuclei of galaxies (around 
10 5 M© - 10 9 M Q ). Although at present there is no decisive astrophysical evidence 
for the existence of compact objects with event horizons, the accurate PN equations 
of motion, as well as the emitted gravitational waveforms, are in principle capable 
of revealing the presence or absence of the event horizon. This is one of the features 
that makes these studies so exciting. 

The inspiralling compact binaries in the last minutes of their orbital evolution 
arc inherently powerful sources of gravitational waves. A crucial practical issue is 
the number of such systems that may populate a typical galaxy during a given 
interval of time, say, 1 year. Since the detectable part of gravitational radiation 
from such sources is a relatively short-lived event, this number is called the event 
rate. If the event rate were too much low, the prospect of observing these powerful 
sources would not be quite realistic. One would need to survey a huge volume of the 
Universe in order to have some confidence that at least several events occur per 1 
year of observations. As a consequence, the sensitivity of the instrument would be 
required to be extremely high in order to guarantee the possibility of detecting the 
sources located in the most remote parts of this volume. 

Unfortunately, the astrophysical event rates are not very certain . 8 ~ 10 Sometimes 
their evaluation differs by orders of magnitude even in the papers of one and the 
same research group. This is especially true with regard to the binary systems 
involving black holes, for which the observational information is scarce. 11 However, 
it appears that there exists growing consensus at least about the neutron star - 
neutron star (NS-NS) coalescences. The currently most quoted rate 3 x 10~ 5 — 1CP 4 
per year, for a galaxy like our own Milky Way, is at the level estimated very early 
on i' 12-14 and advocated for long time by some groups. 3 

It also follows from general arguments based on binary star evolution and their 
numerical simulation that the black hole - black hole (BH-BH) event rate is expected 
to be only an order of magnitude lower than the NS-NS rate. If this is true, the 
larger total mass of the BH-BH systems, in comparison with the NS-NS systems, 
and therefore larger gravitational wave luminosity, can more than compensate their 
lower event rate, when it comes to the analysis of detectability of these systems. For 
an instrument of a given sensitivity, such as LIGO 15 or VIRGO, 16 the probability 
of seeing a BH-BH coalescence turns out to be higher than the probability of see- 
ing a NS-NS coalescence. Specifically, it was estimated that the BH-BH detection 
rate is likely to be a factor ~ 5 higher than the NS-NS detection rate, despite the 
opposite order of these systems event rates. This justifies the expectation that the 
first detected sources will actually be the coalescing stellar mass black holes. 4 The- 
oretical and numerical studies of inspiralling and coalescing black holes are rising 
as a particularly interesting and important field of study. 
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Even the most powerful expected signals will be not more than at the level of 
noise in the presently available (ground-based) instruments. To extract the signal 
from the noise, and to determine the parameters of the radiating system, one needs 
to know in advance, and as accurately as possible, the theoretical signal templates 
for the incoming waves. It is also assumed that the response of the instrument to 
these waveforms is known with the equally high accuracy. The templates are being 
cross-correlated with the noisy output of the detector. Since the signal from inspi- 
ralling binaries is quasi-periodic, the knowledge of its phase is especially important, 
and the templates must remain in phase with the expected true signal as long as 
possible. The calculation of accurate templates from the compact binaries is now a 
matter of great activity. 

Processes involving supermassive black holes (SMBH) 17 provide more scientific 
oportunities, but add new questions and complications. Gravitational radiation from 
a neutron star or a stellar mass black hole inspiralling into a SMBH contains a wealth 
of information. A detailed waveform allows, in principle, to decide whether the cen- 
tral object is a black hole or some even more exotic object. Indeed, the multipole 
moments of a Kerr black hole satisfy unique relationships as functions of its mass 
and spin, whereas the multipole moments of an arbitrary body are not, in general, 
linked by similar relationships. The information about the multipole moments of 
the central object is encoded in the waveform, thus providing an example of strong 
gravity test with the help of gravitational waves. In practice, one could start from 
evaluation of the quadrupolc moment of the central object (which enters as a pa- 
rameter in the 2PN waveform), and check whether it satisfies the conditions for 
a Kerr black hole. However, the complicated trajectories of infalling masses, the 
uncertain radiation reaction force acting on the body, and the inevitable presence 
of accreation disks surrounding real (in contrast to ideal, theoretical) black holes, 
will make the extraction of astrophysical information not so easy. 

Though most of the current activities in relativistic gravity are related to the 
imminent observation of gravitational waves, pulsar astronomy also has the poten- 
tial of further tests of gravitational theories and more advanced insights into the 
Einsteinian gravity. To mention in this respect are further observations of the dou- 
ble pulsar to eventually measure the moment of inertia of one of the pulsars which 
would give information on the equation of state of neutron stars, 18 and the con- 
struction of SKA (square kilometre array) to detect practically all pulsars in our 
galaxy with a chance of discovery of pulsar-black-hole binaries. 19 

2. Theoretical modelling of compact binaries 

One normally makes a convenient, but approximate, separation of the problem 
into several parts: internal, external, and far zone. This division is justified by the 
hierarchy of characteristic length scales of the binary system: size of the bodies, size 
of the orbit, and wavelength of the emitted gravitational radiation. The internal 
part of the problem is concerned with the bodies themselves, including their sizes, 
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shapes, internal structure and tidal gravitational effects induced by one body on 
another. The external problem deals with equations of motion of the centers of 
mass of participating bodies. As we said, in the external problem the bodies are 
often treated as point particles with given masses. In advanced treatments, the 
point particles are endowed also with spins and higher multipole moments. Finally, 
the far zone problem is mostly the calculation of the emitted gravitational waves. 
Typically this is being done in leading order in terms of the distance to the source. 

Obviously, all three parts of the problem are interconnected. In particular, from 
the near-zone equations of motion, or from the balance relationships equating the 
emitted radiation and the changing orbital characteristics of the system, one finds 
the radiation reaction force and corrections to the Kcplcrian parameters of the bi- 
nary, and then makes further corrections to the waveforms. It should be qualified as 
a remarkable success that the problem of theoretical templates has been worked out 
by successive approximations up to the 3.5PN approximation of general relativity, 
corresponding to the accuracy (v/c) 7 , where v is the relative speed of components 
of the binary. 

The external problem (that is, the motion of the centers of mass) has been 
solved at 3.5PN order independently by three groups, with completely equivalent 
results. One group used the Arnowitt-Deser-Misner (ADM) Hamiltonian formal- 
ism of general relativity 20-23 and worked in a corresponding ADM-type coordinate 
system. Another group used a direct PN iteration of the equations of motion in 
harmonic coordinates. 24-26 Both groups used a description by point particles and a 
self-field rcgularization. The end results of these two approaches have been proved 
to be physically equivalent. 23,25 However, both approaches left undetermined one 
dimensionless parameter at the 3PN order. The appearance of this unknown pa- 
rameter was related with the choice of the regularization method used to cure the 
self-field divergencies of point particles. The completion of the equations of motion 
at the 3PN order was made possible thanks to the powerful dimensional self-field 
regularization, which could fix up uniquely the value of the ambiguity parameter in 
both calculations. 27,28 The third approach 29-31 succeeded in obtaining the equiva- 
lent 3PN equations of motion directly, i.e. by using a "surface integral" method in 
which the equations of motion are written in terms of integrals on surfaces surround- 
ing the compact bodies. This method is applicable for extended compact objects in 
the strong-field point particle limit. Finally, the 3.5PN terms, which constitute a 
1PN relative modification of the radiation reaction force (and are relatively easier 
to derive), have been added in Refs. 32-37 

The far zone problem (that is, the radiation field) has also been solved at the 
3.5PN order, i.e. (v/c) 7 beyond the leading order given by the Einstein quadrupole 
formula, using a particular gravitational wave generation formalism combining mul- 
tipolar expansions with the PN approximation. 38 The crucial step is the computa- 
tion of the binary's quadrupole moment at the 3PN order which has been done 
by a combination of Hadamard's self-field regularization dealing with most of the 
terms, 39 and eventually completed, like in the problem of equations of motion, by 
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dimensional regularization able to fix the value of a few remaining ambiguity coeffi- 
cients. 40 The final 3.5PN templates 42 take into account the values of the ambiguiy 
parameters computed in 27 ' 28,40 (for a review, see 5 ). 

However, there are also some difficult issues involved in this program. One should 
point out that the full analysis of the internal problem has not been done at 3PN 
order and is hard to perform. It has been carried out only for spherically-symmetric 
bodies and to the lowest radiative order (v/c) 5 , i.e. to 2.5PN order in the equations 
of motion. Satisfyingly, it was shown that the compactness parameter, character- 
ising the size of the body L in comparison with its Schwarzschild radius r g , can be 
absorbed into the redefinition (or "renormalization" ) of the body's mass. Therefore, 
at this level of accuracy, the equations of motion are valid for compact objects of 
any size and structure, presumably including black holes. This dependence of the 
results only on integral parameters of the compact body (for example, its mass and 
spin) and independence on its actual shape and internal structure (for example, 
rearrangement of layers of matter) is sometimes called the principle of "effacement" 
of the internal structure. 47 Obviously, the effacement principle cannot be true with 
arbitrarily high accuracy. It has to be violated at some sufficiently high order of PN 
approximations. 

Despite the fact that the point-particle approach has reached a great level of rigor 
and sophistication, it cannot be arbitrarily accurate for real physical objects, even if 
it is consistent mathematically for idealized point particles. We have to worry about 
the magnitude of the finite-size effects, typically in the form of phcnomcnological 
parameters of an extended body and ultimately in the form of the finite gravitational 
radius of a black hole. Even though the PN corrections of very high level of accuracy 
will not be required by observers in the near future, further work is needed for 
precise identification of the limits of theoretical consistency of the PN program 
itself. It is expected, though, that the trouble will not show up too early. Indeed, it 
follows from the simple arguments (see e.g. 5 ' 43 ) that the tidal force in a pair of fluid 
bodies is of the order of k(L/R) 5 (GM 2 /R 2 ), where L is linear size of the bodies, 
R - distance between them, and k is a phenomenological parameter characterizing 
the "elasticity" of the bodies. In the limit where L approaches gravitational radius 
r g , the tidal force is a factor k(v/c) 10 smaller than the Newtonian force, i.e. it is 
formally of the 5PN order and, hence, is very small. This expectation is confirmed by 
explicit calculations of finite-size effects for neutron star binaries 41 and is consistent 
with the relativistic equations of motion for black holes derived by matching of the 
perturbed Schwarzschild solutions. 47 

It should be noted, however, that the above division of the problem into the 
internal and external parts certainly breaks down at the merger phase of compact 
objects. This is especially true for merging black holes and subsequent ringdown 
phase of the combined black hole. This area of study requires new techniques and 
approaches (partially developed, see for example 48,49 ). In particular, the merger and 
ringdown of binary black holes have recently been implemented by numerical tech- 
niques 44-46 Although the amount of gravitational radiation emitted at merger and 
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ringdown phases is relatively modest, it may be observable by advanced detectors, 
opening an exciting era of new discoveries. 

3. Post-Newtonian equations of motion 

The introduction to the problem of testing general relativity with gravitational 
radiation from compact binary inspirals was given by C. Van Den Broek and B. 
S. Sathyaprakash (reported by C. Van Den Broeck). The speaker has emphasized 
that some of the PN terms in the measurable gravitational wave phase arise from 
the scattering of gravitational waves off the gravitational field of the source, in the 
vicinity of the binary. These terms are known as gravitational wave tails. Observa- 
tional checks of the presence of such tail terms (using, for example, the data analysis 
method of 50 ) will be tests of the validity of PN approximations and general relativ- 
ity itself. The contribution of C. Van Den Broek and B. S. Sathyaprakash has also 
discussed some other tests of gravitational theories, which are not necessarily based 
on PN approximations. The range of tested theories includes those with massive 
gravitons. We will briefly review these possibilities below. 

S. Kopeikin spoke about the irrelevancy of the internal structure of gravitating 
bodies (i.e. the effacement principle) in the PN approximations of general relativity 
and scalar-tensor theories of gravity. He argued that in general relativity the effacing 
principle is violated by terms proportional to the rotational moments of inertia of the 
fourth order. In the scalar-tensor theories of gravity the violation begins earlier, by 
the terms proportional to the second order rotational moments of inertia. 51 When 
the effacement principle is violated, the equations of motion of extended bodies 
differ from those of point-like particles. Correspondingly, the emitted waveforms 
are also different. In the limit where the size of the body is taken to be close to 
the Schwarzschild radius r g = 2GAl/c 2 , Kopeikin evaluates that the effacement 
principle is violated in the 3PN approximation in the case of scalar-tensor theories, 
and in the 5PN approximation (terms of the order of (w/c) 10 ) in general relativity. 
The latter statement is in agreement with an earlier conclusion by Damour. 47 

The point-particle approach inevitably encounters the necessity of regulariza- 
tion of the fields diverging at the world lines of the particles. There are two some- 
what different methods of dealing with this problem. One is directly applicable 
in the employed harmonic coordinate system 24-26 while the other is based on the 
ADM formalism. 20-23 Both methods ultimately rely on analytical continuation of 
the equations of motion to the (in general, complex) d ^ 3 dimensions. 27,28 In his 
talk, P. Jaranowski spoke about the dimensional rcgularization of the gravitational 
interaction in the ADM formalism. More precisely he regularized the 3PN Hamilto- 
nian of point masses, in which the field degrees of freedom are reduced using the field 
equations. 27 The speaker argued that the dimensional continuation leads to a finite 
and unambiguous Hamiltonian in the limit d — > 3. He showed that three somewhat 
different methods of computation lead to the same expression for the dimensionally 
regularised 3PN Hamiltonian. This increases confidence in the correctness of the 
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3PN Hamiltonian and the associated equations of motion (an alternative method 
which confirmed the result is 29-31 ). 

G. Faye has derived the equations of motion of spinning compact binaries in- 
cluding the spin-orbit (SO) coupling terms 1PN order beyond the leading-order 
effect. 52 For black holes maximally spinning this corresponds to 2.5PN order. The 
result confirms the previous calculation of Ref. 53 The SO effects up to 2.5PN order 
are also computed in the conserved (Noetherian) integrals of motion, namely the 
energy, the total angular momentum, the linear momentum and the center-of-mass 
integral. The spin precession equations at 1PN order beyond the leading term are 
also obtained. The speaker reported then the computation (using the multipolar- 
PN wave generation formalism 38 ) of the SO contributions in the gravitational- wave 
energy flux and the secular evolution of the binary's orbital phase up to 2.5PN or- 
der. 54 It was shown that the 1PN SO effects are in general numerically larger than 
the spin-spin (SS) effects, in terms of the number of gravitational-waves cycles, even 
though they appear at a formally higher PN order. These results provide more accu- 
rate gravitational- wave templates to be used in the data analysis of rapidly rotating 
Kerr black-hole binaries with the ground-based and space-based detectors. 

L. Gergely discussed the corrections to the 2PN equations of motion of struc- 
tureless point-like particles, which arise when the particles are endowed with spins 
and mass quadrupole moments. If the body is a neutron star with a strong magnetic 
field, then the representing "particle" is endowed also with a magnetic type dipolc 
moment. The speaker has demonstrated a generalized Kepler equation 55 where the 
orbital elements include contributions from the spin-spin, mass quadrupole and 
magnetic dipole interactions. The considered effects nicely fit in the standard Ke- 
pler form with the modified parameters of the orbit. 

The orbital phase of inspiralling binaries with the inclusion of the orbit's eccen- 
tricity was discussed in the contribution of M. Vasuth. He presented the results on 
the change of the mean motion parameter for eccentric orbits and the evolution of 
the orbital phase for circular orbits. All linear effects due to spins, mass quadrupole 
and magnetic dipole moments are included in this derivation. The author has re- 
lated the change of the mean motion and orbital frequency with the radiative energy 
losses. This allows one to derive the time-dependent orbital frequency and phase, 
and to calculate the number of relevant gravitational wave cycles. In the 2PN ap- 
proximation this number includes corrections from spin-orbit (SO), spin-spin (SS), 
mass quadrupole - mass dipole and magnetic dipole - magnetic dipole coupling ef- 
fects. A special attention was paid to the presence of a new self-interaction spin 
term. 56 

M. Tessmer has reported on the accurate and computationally efficient deriva- 
tion of waveforms produced by binaries with arbitrary eccentricity and mass ratio, 
and moving in slowly precessing orbits. 57 The orbital motion is restricted to be 
1PN accurate, and only the quadrupole contribution to the waves' polarization 
amplitudes is taken into account. The central point of the derivation is a special 
numerical method for solving the participating generalized Kepler equation. The 
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speaker has discussed the relevance of the derived waveforms for observations with 
the space-based interferometer LISA. 61 

G. Schafer has reported on recent results from his research group additionally 
to Tessmer's contribution, which were the phasing of gravitational waves from in- 
spiralling eccentric binaries at the 1PN radiation-reaction order (corresponding to 
3.5PN order in the equations of motion) 58 and the gravitational recoil during binary 
black hole coalescence using the effectivc-onc-body approach. 59 

L. Blanchet discussed the problem of the gravitational recoil of black-hole bi- 
naries using PN techniques. He reported on the recent calculation of the recoil of 
non-spinning black holes in the inspiralling phase at 2PN order. 60 The author esti- 
mated the kick velocity accumulated during the plunge from the innermost stable 
circular orbit up to the horizon by integrating the momentum flux along a plunge 
geodesic of the Schwarzschild metric. The contribution to the total recoil due to the 
subsequent ringdown phase was neglected. 



4. Supermassive black holes and strong field tests of gravity 

In addition to possible SMBH mergers, a promising source for LISA, called the 
extreme mass ratio inspiral (EMRI), is the infall of a compact object into a SMBH. 
The infalling body can be a neutron star or a stellar mass black hole. The detailed 
study of gravitational waves emitted by an EMRI event can, in principle, allow 
one to build a "map" of the gravitational field in a parscc radius of the galactic 
nucleus. It is presumed that the SMBH candidate in the nucleus is a Kerr black hole. 
However, the modelling of gravitational wave from EMRIs is a very challenging task 
as the orbits often exhibit complicated behaviour, and the EMRI parameter space 
is huge. In this situation, the derivation of true waveforms is practically impossible 
by analytical techniques and is very expensive by numerical methods. 

One approach to the problem is the construction of a family of simpler, approx- 
imate waveforms, called "kludge" waveforms, which nevertheless capture the main 
features of the true signals. 62 S. Babak has discussed the progress in compiling a 
bank of detection templates which are numerical kludge (NK) waveforms. They 
can be generated quickly and cheaply, and they are good enough to be used as 
a first pass test for the parameter estimation. It is reported that satisfactory NK 
waveforms are available for infalling masses in the range of up to 5-6 M Q . 

The derivation of sufficiently accurate waveforms always requires the proper 
taking into account of the radiation reaction force acting on the body and changing 
its trajectory. This problem is especially acute for EMRI sources where the inspiral 
phase lasts for long time and orbits are complicated. In SMBH studies, one normally 
considers the technique of small perturbations of the background space-time. Due 
to the emission of gravitational waves, the world line of a compact object deviates 
in a calculable way from the nominal geodesic line of the background (possibly, a 
Kerr black hole) solution. 63 

In the point-particle approximation, the metric perturbations are divergent along 
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the particle's world line. This requires some sort of field regularization. Although 
many technical problems concerning the regularization are already solved, 64 some 
conceptual issues remain. Y. Mino has reported on the existing problems for the 
Kerr-background metric related to the choice of gauge, the restrictive assumption 
of linearity (weakness) of perturbations, and finally, the extraction of gravitational 
waveforms. One way to resolving these difficulties is a careful identification of the 
radiation-reaction part in the self-force expression. The author has also discussed 
the guides provided by the adiabatic approximation which assumes that the orbital 
parameters evolve slowly. 

We normally take it for granted that compact binary systems are "clean", so 
that the relativistic gravity is essentially the only participating interaction up to the 
very late stages of evolution near the final merger. However, realistic astrophysical 
compact objects are likely to be surrounded by accretion disks which will complicate 
the motion of close companions and may intervene at the level of corrections larger 
than the magnitude of relativistic PN effects. This seems to be especially plausible 
to happen in the case of EMRI sources, and this was the subject of the contribution 
by P. Basu. It is indeed expected that the late orbits of the infalling object will be 
taking place within the accretion disk of the SMBH. The infalling NS or BH will 
itself accrete some matter from the disk. 

Typically, the specific angular momentum of the accreted matter is lower than 
the Keplerian angular momentum of the infalling body. Therefore, the total angular 
momentum of the infalling object will decrease more rapidly than expected, leading 
to a faster infall into the SMBH. For some accretion disks the situation can be 
opposite, which would lead to the slower than expected infall of the companion. 
Basu has outlined the gravitational field and fluid dynamics equations that should 
be solved simultaneously. The gravitational field of the central Kerr black hole is 
modelled by an effective Newtonian potential. 65 For disks with certain parameters, 
the change of the companion's angular momentum can be greater than the largest 
losses expected due to gravitational radiation. Qualitatively, this situation is similar 
to what can happen in neutron star systems when interactions other than gravity 
are present. 4 

5. New research directions 

The overwhelming majority of current studies in gravitational physics are based 
on the geometrical formulation of Einstein's equations. For example, in PN expan- 
sions describing an isolated binary system (usually, in asymptotically Lorentzian 
harmonic coordinates), one normally treats the gravitational functions /i Miy (i,x) as 
frame-dependent pieces of curved space-time metric g^v^, x), rather than compo- 
nents of a genuine tensorial gravitational field defined in a flat space-time with 
Minkowski metric. 

The thrifty geometrical picture of general relativity, which combines gravity, ge- 
ometry and a choice of coordinates in a single mathematical object - the curved 
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space-time metric tensor g^(x a )- is accompanied by some peculiar features. One 
may recall ambiguities in the description of the gravitational field energy density 
and gravitational- wave fluxes, constant mixing of external (coordinate) and inter- 
nal ("gauge") transformations, and, in general, certain disjointment of geometrical 
gravity from other field theories. 

It is known for long time that the geometrical picture of general relativity is by 
no means compulsory or necessary. It is argued 67 that the geometrical Einstein's 
gravity is fully equivalent mathematically and physically (including cosmology) to 
a field theory in a flat space-time, i.e. in a space-time with zero-curvature metric 
tensor j flu (x a ). In flat space-time, one can always choose global Lorentzian coordi- 
nates, so that the metric tensor j^ix 01 ) takes on the familiar form of the Minkowski 
metric rf ,v . The field-theoretical formulation of general relativity possesses all the 
necessary and strictly defined structures: covariant gravitational Lagrangian con- 
taining gravitational field h^ v (x a ) and its first derivatives, second-order differential 
non-linear field equations, gravitational energy-momentum tensor free of second 
and higher order derivatives of the field, universal coupling of the gravitational field 
to other physical fields (which allows one to combine W{x a ) and ^^(x a ) into a 
single object g flv (x a ) and to reinterpret the theory as a geometrical theory, where 
the curved space-time metric tensor g^ u (x a ) satisfies geometrical Einstein equa- 
tions), conservation laws, well defined and physically distinct coordinate and gauge 
freedoms. 

One does not expect that any new observational conclusions will arise from a 
reformulation of one and the same fundamental theory, but a new angle of view 
helps one to see problems in a different light and answer questions which otherwise 
could not be even properly formulated. In particular, the field-theoretical approach 
to gravity opens the door for natural modifications of general relativity including 
the concept of massive gravitons. 67 It is not surprising that what we are doing in 
practice, "by hands" , in PN expansions is similar to traditional field-theoretical per- 
turbative calculations, independently of whether we adhere to geometrical ideology 
or not. 

The relevance of field-theoretical techniques for PN calculations in binary sys- 
tems was emphasized by /. Rothstein and R. Porto. It was argued that the replace- 
ment of detailed internal structure of the gravitating bodies by a set of phenomcno- 
logical parameters associated with a point-like particle is similar to what is routinely 
being done in effective field theories (EFT) where one is interested in a simple de- 
scription of the influence of short-scale physics on large-scale (low energy) dynamics. 
The well developed methods of EFT allow one to introduce simple power-counting 
arguments (evaluation of the order of magnitude of various contributing terms) and 
employ traditional treatments of divergencies and renormalization procedures. 68 ' 69 
The formalism also allows one to calculate absorptive effects for an arbitrary object 
in terms of the graviton absorptive cross-section. 68 (In the context of gravity, the 
authors are using a somewhat confusing name of "Non-Relativistic General Rela- 
tivity".) For application to PN computations of compact binaries, these methods 
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should be applied twice. First, when one derives the equations of motion for point 
particles, and second, when one calculates the gravitational waveforms and replaces 
the entire binary system by a single particle with certain multipole moments. The 
back-reaction effects within the EFT approch should also be properly worked out. 

/. Rothstein has stressed that the EFT calculations are essentially being done at 
the level of the Lagrangian and the action. He presented the action for a binary with 
slowly moving components. He argued that the divergencies arising at the (v/c) 6 
level, that is, in the 3PN approximation, are not physical as the corresponding terms 
in the action can be removed by field redefinitions. This conclusion is consistent with 
what was stated on the grounds of calculations in the more usual PN approach. 24,28 
Interestingly, the same line of EFT arguments has led the speaker to the conclusion 
that first terms which cannot be removed by field redefinitions, and therefore provide 
a source of violation of the "effacement" principle, are of the order of (v/c) 10 , i.e. 
they appear in the 5PN approximation. Again, this conclusion is consistent with 
other arguments, thus increasing our confidence in the internal workings of the 
entire PN scheme. 

R. Porto has extended the EFT approach to include spin dynamics. This is 
achieved by adding rotational degrees of freedom to the world-line action of the 
point-like particles. The issue of different choices for the spin supplementary condi- 
tions was clarified. It was shown that these conditions are equivalent, in the sense 
of the final results, at least at the 1PN level. In the areas where the EFT approach 
overlaps with previous studies, 70-73 the final conclusions are in agreement with each 
other. The author has also reported new results, derived by EFT techniques, on the 
corrections to the spin-spin (SS) potential in the 3PN approximation. 74 It is argued 
that these corrections are easier to handle in the EFT approach, and, in general, that 
this approach is a powerful tool to treat in a systematic manner the higher-order 
PN effects. 

Spin effects encapsulated in a prescribed Lagrangian were also considered in 
the contribution of M. Vasuth. It is assumed that the motion of the binary system 
is described by the Lense-Thirring Lagrangian which takes into account the spin 
vector of the central body and treats another body as a spinless test particle. The 
author considers the 1.5PN accurate motion of the test particle and concentrates 
on terms linear in the spin of the central body. The radial and angular dynamics 
of the system are treated with the help of more general results. 75 The outcome of 
calculations are both polarization components of the emitted gravitational waves, 
including the contributions linear in the spin. It is confirmed that the waveforms 
are in agreement with previous calculations. 71 

The geometrical general relativity requires one to be careful with such things 
as choice of coordinates, identification of the "gauge-independent" degrees of free- 
dom of the gravitational field, clock synchronization, equivalence principle, and, in 
general, physical interpretation of the participating quantities. 76 L. Lusanna has 
discussed some of these issues in his contribution. Although these issues refer in 
general to the full theory of geometrical gravitation, they also manifest themselves 
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in approximations, such as PN gravitational-wave studies. 

Returning to the review of tests of general relativity, it is important to remem- 
ber that interesting tests are not necessarily associated with the strong field regime. 
C. Van Den Broek and B. S. Sathyaprakash have reminded us of the importance 
of some signatures in the regime of weak plane gravitational waves. According to 
general relativity, the response of an interferometer to the incoming plane wave 
contains not only the usual "electric" component but also the (typically, smaller) 
"magnetic" component. The identification of the "magnetic" component in the out- 
put data is needed for proper extraction of the radiating system's parameters, but 
also as a test of the "magnetic" prediction of general relativity. 66 

New possibilities arise in alternative theories of gravity, such as theories with 
massive gravitons. In particular, the detection of a "scalar" polarization state of 
gravitational waves, expected in such theories, would revolutionize our views on the 
nature of gravity and could possibly provide a gravitational explanation to some 
presently existing cosmological puzzles. 

6. Conclusions 

The current advances in observational facilities (LIGO and VIRGO on Earth, LISA 
in Space) have stimulated further deep insights in such problems as gravitational- 
wave physics in general and sources of gravitational waves in particular, relativistic 
celestial mechanics, tests of general relativity and alternative theories. The con- 
tributions to the PPN1 and PPN2 parallel sessions have demonstrated the depth 
and rigor of the continuing research in these areas. Although several of the recently 
derived results need to be cross-checked and placed in a unique common context, 
it is clear that there is no major conceptual or technical difficulties in this field. In 
particular the status of PN computations of equations of motion and gravitational 
radiation is quite satisfactory. In the coming years we will probably witness a major 
progress in relating and connecting the analytical PN calculations with successful 
numerical computations. However, it appears that there is already enough theoreti- 
cal clarity and completeness at least for the level of accuracy of existing gravitational 
wave observations. Hopefully, further theoretical work will proceed hand in hand 
with successful experiments. 
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